#include "pdeDoseEngine.h"

using namespace pde;

void main() {
    std::string path1("f:\\proton_dose_models\\test1\\");
    auto nozzle = Nozzle::New();

    nozzle->from_file(path1 + "nozzle.bson");

    auto dose_model = nozzle->get_dose_models()[0];
    std::vector<double> z_list = {-200.0, -100.0, 0.0, 100.0, 200.0};

    auto sigma_list = calculate_air_profile_sigma(dose_model, z_list);
}